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The ALEPH detector at LEP has been used to measure the momentum spectrum and charge ratio 
of vertical cosmic ray muons underground. The sea-level cosmic ray muon spectrum for momenta 
up to 2.5TeV/c has been obtained by correcting for the overburden of 320 meter water equivalent 
(mwe). The results are compared with Monte Carlo models for air shower development in the 
atmosphere. From the analysis of the spectrum the total flux and the spectral index of the cosmic 
ray primaries is inferred. The charge ratio suggests a dominantly light composition of cosmic ray 
primaries with energies up to 10^^ eV. 



The sea level cosmic ray muon spectrum resulting from 
the interaction of primary cosmic rays with the atmo- 
sphere probes multi-particle production in hadronic in- 
teractions up to and beyond LHC energies. In addition, 
the muon charge ratio is sensitive to the fraction of heavy 
nuclei in the primary cosmic ray beam. 

In the past many measurements were based on sim- 
ple detection devices at the earth's surface, see e.g. pj 
and references therein. Additional information can be 
obtained from the highly sophisticated detectors used in 
particle physics |2H5]. The fact that those detectors are 
placed under ground offers the additional advantage of a 
natural muon filter and a momentum cutoff which em- 
phasizes the high energy part of the spectrum. Here we 
present a measurement of the vertical muon spectrum 
and charge ratio performed with the ALEPH detector at 
LEP, which was located at a depth of 320 mwe. 

The ALEPH detector is described in detail in . The 
study presented here uses mainly the large time projec- 
tion chamber (TPC), the main tracking device of the 
ALEPH detector. The data were taken in 1999 during 
one week of dedicated running with a cosmic trigger us- 
ing information from the hadron calorimeter (HCAL). 
Background from interactions in the iron of the HCAL 
and electromagnetic calorimeter (ECAL) can be reliably 
identified, have been shown to be well described by the 
detector simulation |7] and are taken into account in the 
Monte Carlo corrections. 

The coordinate system used in the analysis is the stan- 
dard ALEPH coordinate system. It is a right-handed 
system with the cc-direction pointing towards the center 



of LEP, the y-direction being vertically upwards and z 
along the beam-direction. The center of the detector is 
'At X = y = z = 0. Muon tracks are selected by the 
requirement of having at least 27 TPC hits out of a max- 
imum of 42, a reconstructed momentum p > lOGeV/c 
and a polar angle with respect to the vertical 9 < 20°. 
The cut on the number of TPC hits guarantees precise 
momentum reconstruction, the minimum momentum cut 
ensures that the muon has enough energy to pass through 
the entire barrel of the calorimeter and thus full trigger 
efficiency. Both cuts reject low energy background from 
secondary interactions. The quality of the reconstruction 
is demonstrated by the fact that such details of the geo- 
metrical structure as the access shaft of the ALEPH pit 
and the flux shadowing by the Jura mountains are visible 
in the arrival directions of the incident muons. For the 
final analysis additional fiducial cuts were applied. Az- 
imuthal angles (measured around the vertical) of ±10° 
with respect to the direction of the underground access 
shaft were cut out to mask the enhanced muon flux from 
this direction. In addition, the muon tracks were required 
to pass through the central plane y = in the range 
|a;o| < 50 cm, |zo| > 50 cm and l^ol < 170 cm. These cuts 
make sure that the selected tracks are fully contained in 
one half of the TPC and pass through the full diameter of 
the chamber, optimizing momentum resolution and min- 
imizing systematic effects due to boundary effects and 
the inversion of the drift field at the HV-membrane at 
z = 0. The momentum resolution dp/p of the TPC is 
shown in Fig. [T] as a function of p for three values of the 
fiducial cut on xq. One clearly sees how the performance 
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improves towards the central region. The nominal cut 
value \xo\ < 50cm is chosen as a compromise between 
good resolution and still high statistics. Here the maxi- 
mum measurable momentum is around 2TeV/c. Extrap- 
olated to the surface, the accessible momentum range is 
up to 2.5TeV/c. 



momentum resolution of the TPC 



1 



10' 

















! 






























































































■ */ 
































































1 




.*/ 




























: 


y 


































*, 








































































|x |<70 cm 
|x^|<50 cm 
|xj<30 cm 
































































































































































Kf l(f 
p[GeV/c] 



FIG. 1. Momentum resolution as a function of the momentum 
for three cuts on the fiducial area of the TPC. 

The raw distributions for spectrum and charge ratio 
were corrected using momentum dependent correction 
factors obtained from a full simulation of single muons 
passing through the ALEPH detector from above. For 
this simulation, at a depth of 320 mwe, an energy spec- 
trum proportional to was assumed, with 7 — —3. 
The correction factors, parameterized by 4th order poly- 
nomials as a function of Igp, include effects due to elec- 
tromagnetic interactions of muons in the detector ma- 
terial, multiple scattering in the calorimeters and finite 
momentum resolution of the TPC. For momenta from 10 
to 40 GeV/c the correction factors for the spectrum vary 
from values around 2 to 1.2. Towards higher momenta 
they decrease slowly, falling below unity around 1 TeV/c. 
For the charge ratio no trend is observed. With aver- 
age corrections below 1 %, the analysis of the charge ra- 
tio was performed without explicit detector corrections. 
In the end a global 1 % systematic error, correlated be- 
tween all bins, was assigned. The nominal results for 
the muon momentum spectrum and charge ratio at the 
surface were then determined by propagating the mea- 
sured momenta to the surface, using a parameterization 
of the form dE/dx — a + bE with energy dependent 
coefficients a and b [S]. The continuous energy loss is 
approximated by performing the propagation in five dis- 
crete steps, which has been checked to be within 1% of 
the continuous treatment. An overburden of 306 mwe 
for (j) < 180 deg (towards lake Geneva) and 326 mwe for 
(p > 180 deg (towards Jura mountains) was assumed, 
which takes into account the surface topography close 
to the ALEPH experiment. After propagation to the 



surface the spectrum measured up to zenith angles of 
20 degrees was corrected to the purely vertical spectrum 
by a momentum dependent factor I{p,0)/I{p,9) using 
the parameterization of the muon flux in reference [S]- 
The correction is in the range 1-3%. The absolute flux 
finally is obtained from 

^corr 



Here N^°'^^ is the number of muons measured after all 
corrections in the momentum interval [p,p + Ap], S is the 
fiducial area in the TPC and ft the solid angle accepted 
for the measurement. The uptime T is estimated from 
an independent measurement [lU] of the CosmoALEPH 
HCAL trigger rate and the number of muons per unit 
area observed during the data taking as T = (3.52 ± 
0.32) • 10^ s where the error is dominated by systematics 
estimated from spatial inhomogeneities of the observed 
muon flux through the central plane of the TPC. Details 
of the analysis are described in references [H [TT]. The 
total normalization uncertainty is 9.1%. 

For the determination of the systematic uncertainties 
of the results the analysis was repeated 10** times, ran- 
domly varying global settings of the analysis and fluc- 
tuating the individual measurements within their uncer- 
tainties. With D2°"^ the result from the nominal analy- 
sis, a covariance matrix describing the systematic uncer- 
tainties and correlations between all bins was constructed 
from those pseudo-experiments by 

= ((A - OmiDj - D^n) ■ (2) 

Individual tracks were subjected to smearing of 1/p 
according to the estimated error from the track fit, plus 
direction and position smearing due to multiple scatter- 
ing in the overburden and the finite resolution of the TPC 
as parameterized from Monte Carlo simulations. Global 
variations between different pseudo-experiments which 
have been considered are the spectral index in the calcu- 
lation of the detector corrections, 7 G {—2.7, —3.0, —3.3}, 
the use of an additional quality cut dp/p < 1 in the track 
selection, omitting the zenith angle correction, or, in or- 
der to probe the sensitivity to the simple energy loss 
model, doing the propagation through the overburden 
in a single step. Continuous global parameters entering 
the analysis have been randomly varied between pseudo- 
experiments by assuming a uniform distribution within 
their assumed uncertainties. The overburden was var- 
ied by ±3.4 mwe on the lake side and ±13.4 mwe on the 
Jura side, where the difference accounts for uncertainties 
in the exact effects of the surface topology. Momentum 
dependent biases on position, direction and momentum 
measurements in the TPC are assumed to be taken into 
account by the Monte Carlo simulation. Nevertheless, 
in order to probe the sensitivity to biases seen in the 
simulation, these were parameterized and added with a 
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p [GeV/c] 
112 
141 
178 
224 
282 
355 
447 
562 
708 
891 
1122 
1413 
1778 
2239 



spectrum [(GeV/c 
(1.959 ± 0.018 ±0. 
(9.004 ± 0.108 ±0. 
(4.519 ± 0.063 ±0. 
(2.202 ± 0.040 ±0. 
(1.072 ± 0.025 ±0. 
(5.068 ± 0.152 ±0. 
(2.373 ± 0.090 ±0. 
(1.204 ± 0.058 ±0. 
(5.745 ± 0.034 ±0. 
(2.745 ± 0.214 ±0. 
(1.209 ± 0.126 ±0. 
(6.741 ± 0.829 ±0. 
(2.217 ± 0.419 ±0. 
(1.536 ifc 0.307 ±0. 
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charge ratio 
1.252 ±0.024 ±0.046 
1.293 ±0.030 ±0.025 
1.259 ±0.036 ±0.017 
1.271 ±0.046 ±0.028 
1.239 ±0.057 ±0.033 
1.348 ±0.080 ±0.048 
1.541 ± 0.121 ± 0.085 
1.373 ±0.133 ±0.101 
1.243 ±0.153 ±0.127 
1.547 ±0.248 ±0.168 
1.785 ±0.388 ±0.369 
1.361 ±0.339 ±0.368 
0.648 ±0.251 ±0.805 
1.495 ±0.610 ±0.698 



TABLE I. Vertical muon spectrum and charge ratio at sea 
level for momenta between lOOGeV/c and 2.5TeV/c. The 
first column gives the momentum at the bin center, the second 
and third columns are the vertical flux and the charge ratio, 
respectively. The first error is statistical, the second one the 
systematic uncertainty of the results. In addition there is a 
global normalization uncertainty of 9.1% for the spectrum. 
The correlation matrix of the systematic errors is given in a 
supporting document to this analysis 



scale factor randomly chosen from a uniform distribution 
in the range [—1, 1], i.e. covering the range between zero 
and twice the estimated bias. In the same way also the 
impact of the uncertainty of the global momentum scale 
cr(l/p) = 3.7 • 10-VGeV/c [12] was studied. 

The CosmoALEPH results for the vertical muon spec- 
trum and the charge ratio are listed in table |l] and dis- 
played together with predictions from different interac- 
tions models in Fig. [2] The spectrum agrees with the 
global parameterization given in reference |T]. Taking 
the full covariance matrix of the result into account, one 
finds a x^-confidence level p = 18.6 %. At large momenta 
the CosmoALEPH measurements are also in good agree- 
ment with the results by the MARS collaboration [13], 
which was historically the first to measure momenta be- 
yond ITeV/c. 

Within errors the measurements of the charge ratio are 
consistent with being independent of the momentum, but 
also with an increase at large energies as observed e.g. 
in [3H5]. Fitting a constant value gives a x^-confidence 
level p = 0.85 and the charge ratio = 1.289 ± 0.022, 
which is in good agreement with the result from [J based 
on momenta up to 316GeV/c. The CosmoALEPH mea- 
surement reaches up to 2.5TeV/c, albeit with large un- 
certainties for momenta above ITeV/c. 

The CosmoALEPH measurements were compared with 
Monte Carlo simulations using the models DPI^JET 14J, 
VENUS 15J, NEXUS [16j, EP0SL17 , QGSJET^IS , QGSJETII^lQj 
and SYBILL HSOj from the CORSIKA program suite, ver- 
sion 6.375 |21j . Studies of different parameterizations 
of the primary spectrum ^22j showed that the data are 
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FIG. 2. Measured momentum spectrum and charge ratio com- 
pared to best fit Monte Carlo predictions. At the bottom the 
error bands of the Monte-Carlo curves are indicated. 



reasonably well described by the Constant Mass Com- 
position model (CMC), both in shape and absolute nor- 
malization. The momentum range covered by the Cos- 
moALEPH data is found to be insensitive to the exis- 
tence of the knee in the primary spectrum. The observed 
charge ratio depends mainly on the average charge per 
nucleon of the primary particle initiating a shower, i.e. it 
is sensitive to the fraction of heavy particles in the pri- 
mary spectrum, but docs not provide information about 
the chemical composition of those nuclei. The momen- 
tum spectrum of cosmic ray muons is sensitive to the 
total fiux and the spectral index of the primary spec- 
trum. 

To determine the fraction of heavy nuclei jh , total flux 
and spectral index 7, Monte Carlo predictions using sim- 
ple power laws 



dn 
dEh 



and 



dn 
dK, 



(3) 



for the heavy particles and protons were fitted to the 
data. For heavy nuclei the cocktail of the CMC-model 
with about 30% He-, 20% N-, 30% Mg- and 20%Fe-nuclei 
was assumed. Technically, in the fit the existing simula- 
tions were reweighted to conform to eq.|3] The parame- 
ters fh and 7 were scanned in the range 2.3 < 7 < 2.9 
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FIG. 3. x^"Confidence level contours as a function of 7 and 
fh- The plot shows the combined confidence level contours 
from the models with a good global fit. 



ray muons at the surface in the momentum range from 
lOOGeV/c to 2.5TeV/c. The charge ratio is sensitive to 
the fraction of heavy primaries. Based on simulations 
with different interaction models, the fraction of heavy 
nuclei is inferred to he fh < 41%. The preferred spectral 
index for the primary spectrum up to energies of 10^^ eV 
is found to be in the range from 2.47 < 7 < 2.59, and the 
flux of cosmic ray primaries with energy above 10^ GeV 
is determined as $2 = 3.74 ± 0.61 m^^s^^sr^^. 

The authors gratefully acknowledge the help of the 
ALEPH collaboration, in particular Markus Frank and 
Beat Jost, as well as Alois Putzer, Bertram Rensch, 
and Thomas Ziegler in doing the measurements. The 
analysis of the CosmoALEPH experiment has been sup- 
ported by the Deutsche Forschungsgemeinschaft under 
grant DFG/Gr/1796/1-3. 



and < //i < 1. At each point the normalization K of 
the spectrum was adjusted such that the model repro- 
duces the best measured point of the momentum spec- 
trum. Then the total for the remaining 27 degrees of 
freedom from spectrum and charge ratio was calculated. 
The two distributions were assumed to be independent, 
but the full correlation matrix within each was properly 
taken into account. The best fit results are shown in 
Fig. [2] 

One finds that all models except SYBILL achieve a good 
description of the shape of the muon spectrum. Good 
global fits are obtained with the DPMJET, EPOS, NEXUS 
and VENUS models. The others fail to reproduce the ob- 
served charge ratio. The four models with a good global 
fit favor a predominantly light composition of cosmic pri- 
maries. For a quantitative conclusion the results from 
those models were combined such that for every point in 
the (7, //i) plane the smallest of the four x^-confidence 
levels was taken as the combined one. This procedure 
yields regions in the parameter plane, where all models 
under consideration are consistent with the data with at 
least the specified confidence level. The result is shown in 
Fig. m Taking the region with common confidence level 
p > 0.1 finally yields 2.47 < 7 < 2.59 and fh < 0.41. 

The fiux normalization obtained from the fit is strongly 
correlated to fh and 7. Expressed through the total 
flux $2 of cosmic ray primaries with energies above 
10^ GeV and considering only the region with common 
confidence level p > 0.1, one finds $2 = 3.68 ± 0.51 ± 
0.33m~^s~^sr~-'^. The central value is the average of the 
results from the models with a good fit to the data, the 
first error the RMS spread of those results and the sec- 
ond one the 9.1% normalization uncertainty of the Cos- 
moALEPH measurement. The statistical errors of the 
result are negligible. 

In summary, we have presented a precision measure- 
ment of the spectrum and charge ratio of vertical cosmic 
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